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It is quite plausible that the mass of the dark matter particle increases significantly after 
its freeze-out, due to a scalar field rolling to large values. We describe a realization of this 
scenario in the context of thermal inflation which naturally gives a cold dark matter particle 
with the correct cosmological abundance and a mass around I0 10 GeV, evading the conven- 
tional upper bound of 10 5 GeV. We also discuss another realization which could produce 
a cosmologically interesting abundance of near Planck mass, possibly electromagnetically 
charged, particles. The detection and observational consequences of superheavy cold dark 
matter or WIMPZILLAs are briefly examined. 

I. INTRODUCTION 

A 'model-independent' bound of about 10 5 GeV 0] has often been invoked (see and references therein) 
as the largest possible mass the dark matter particle can have. The derivation of this bound uses the unitarity 
bound on the annihilation cross-section, and makes one crucial assumption: thermal equilibrium in the early 
universe. The unitarity bound on the annihilation cross-section tells us 

M < ± (i) 

where M is the mass of the dark matter particle, and Jxr^v) is the thermally averaged annihilation cross- 
section that appears in the relevant Boltzman equation 0] (i.e. annihilation per unit time is given by n(o~Av), 
where n is the proper number density of the dark matter particle). The assumption of thermal equilibrium, 
on the other hand, tells us that the freeze-out abundance is given by the thermal distribution 

n~(M7»iexp(-|^) (2) 



where Tf is the freeze-out temperature. Q The exponential suppression implies that Tf cannot be too much 
smaller than M; hence, combining this with the unitarity bound, it can be shown that Q| 
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: This assumes a cold relic. For a hot relic that freezes out at Tj > M, the number density will be higher leading to 
a stronger bound on its mass M. 
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^ 01 (to4v) < 3 » 

where f2 is the ratio of the mass density in the dark matter particle to the critical density today, and h is the 
Hubble constant in units of lOOkms Mpc . This bound naturally has important implications for dark 
matter searches. 

Recently, it has been shown that this bound can be evaded by violating the assumption of thermal 
equilibrium in the early universe, for instance by producing the dark matter particle at the end of inflation via 
preheating/reheating or gravitational particle-creation The masses required for a present abundance 

of fl ~ 1 are within a few orders of magnitude of the Hubble parameter at the end of inflation, i.e. ~ 
10 12 — 10 16 GeV. For this reason, they have been called superheavy dark matter or WIMPZILLAs. 

Here, we propose a different production mechanism that can also evade the 10 5 GeV upper bound, and 
naturally achieves Q ~ 1. It makes use of a late period of inflation called thermal inflation ||-^| that occurs 
at an energy scale of around 10 6 GeV, with a Hubble parameter of the order of 1 keV, which is to be compared 
with the GUT-scale ordinary inflation having V 1 / 4 - 10 16 GeV and H ~ 10 13 GeV used in |J|. 

Thermal inflation provides a natural solution to the Polonyi or moduli problem |J that generically arises 
in string theory (see also jlO]]). It occurs when a 'flaton', a scalar field with a small mass and a large 
vacuum expectation value, is held at the origin by its finite temperature effective potential. One gets a few 
e-folds of inflation because the flaton's potential energy dominates the thermal energy density well before the 
temperature drops below the critical tempe ratur e for the flaton to start rolling away from the origin. The 



prototypical flaton potential is described in § II A . Thermal inflation occurs at a very low energy scale which 
is the reason it can successfully dilute the potentially harmful moduli produced after ordinary inflation. Note 
that it will also dilute any superheavy dark matter produced at the end of ordinary inflation. 

Our idea for dark matter production works roughly as follows. A particle ip and its antiparticle ip, which 
carry some conserved charge to make them stable, are coupled to the flaton. They are initially massless 
during the thermal inflation when the flaton is held at the origin by the finite temperature effects of ip and 
ip. After the temperature of the universe drops below the critical temperature, the flaton begins to fast-roll 
down its potential. The ip particle quickly gains mass in the process, which reduces its annihilation cross- 
section, and its abundance quickly freezes out. After that, the flaton continues to roll until it reaches the 
true vacuum, acquiring a large expectation value and giving ip a large mass.^ The parameters for thermal 
inflation give a mass for ip of around 10 10 GeV, and work out naturally to give an abundance of ~ 1. 
The details are explained in §H 

This production mechanism evades the conventional upper bound of 10 5 GeV by giving the particle ip a 
larger annihilation cross-section at freeze-out than what one would expect based on its final mass, and by 
entropy production after the freeze-out. 

Since thermal inflation provides a natural solution to the moduli problem, and since the prediction of 
~ 1 follows rather naturally from its parameters (assuming ip is stable), the possibility of a sig nific ant 



fraction of the universe being composed of superheavy dark matter should be taken seriously. In § III , we 



2 A related mechanism has also been considered in JTTJ . 
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discuss the observational consequences and detectability of such dark matter, which could be completely 
inert, weakly-interacting, or even electromagnetically or strongly-charged. Finally, we conclude in §IV with 
discussions of other plausible realizations of the mechanism outlined above, which include the possibility of 
producing near Planck mass relics that could perhaps be electromagnetically charged. 



II. SUPERHEAVY DARK MATTER FROM THERMAL INFLATION 



A. Particle physics 

A simple model of thermal inflation is provided by the superpotential [ fL2| 



W(<t>) 



4Mpi 



(4) 



where <fi is the Baton, and Mp\ ~ 2.4 x 10 GeV. This form for the superpotential can be guaranteed by a 
discrete gauge symmetry because the superpotential is a holomorphic function of cj>, i.e. does not depend 
on ^'s complex conjugate (p' . The supersymmetric part of the scalar potential is then given by 



V susy (^) = 



dW 



d<p 



(5) 



In addition, one requires (/>'s soft supersymmetry-breaking mass-squared to be negative. The scalar potential 
is then 



V(^)=V ti -ml\<f>\ 2 +(^ 



+ c.c. 



(6) 



where Vtj and A are other soft supersymmetry-breaking parameters. One would expect \A\ < m^. The scale 
of the soft supersymmetry-breaking parameters of the Minimal Supersymmetric Standard Model (MSSM) 
is expected to be around 100 GeV to ITeV. This potential has four degenerate minima ^ with \<j>\ = vac 
where 



6 2 

vac 




(7) 



3 The domain walls associated with this degeneracy are most likely harmless, either because the vacua are identified 
because of a discrete gauge symmetry, or because higher order terms, that would generically be present in the absence 
of a gauge symmetry, break the degeneracy. 
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For ?7i0 = 100 GeV, \\#\ — 1, and neglecting A, this gives (f> vac — 10 10 GeV. Requiring zero cosmological 
constant at the minima gives 
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(8) 



Thermal inflation starts when the energy density of the universe starts to be dominated by V tl: and ends 
a few e-folds later when the temperature drops below that required to hold (p at <p — 0. <p then rapidly 
rolls towards, and oscillates about, the minima of its potential. It eventually decays, leaving a radiation 
dominated universe, at a temperature Idee- We require 

T dcc >10MeV (9) 

to avoid interfering with nucleosynthesis. 

In order for <p to be held at (p = by thermal effects during the thermal inflation, <p must have unsuppressed 
interactions with other fields in the thermal bath. We will assume these interactions include a coupling of 
the form ^ 

W = X^cpipip (10) 

with | A,/, | ~ 1. After thermal inflation, ip and ip will acquire masses 

M=|A^|0 vac (11) 

from this coupling. For |A^| = 1 and vac = 10 10 GeV, this gives M = 10 10 GeV. 

In order to satisfy the decay constraint, Eq. (0), we require (p to have some, possibly indirect, couplings 
to the MSSM, beyond the ever present gravitational couplings. We will assume this is achieved by having 
ip and ip charged under SU(3)xSU(2)xU(l). Other possibilities were considered in Ref. 0. In order not to 
interfere with gauge coupling unification, ip and ip should form complete representations of SU(5). These 
couplings give a decay rate jl3| 

r,~ 3 x io- 9 ( (12) 

where is the number of internal degrees of freedom of ip and ip. For example, if ip = 5 and ip = 5 then 
=40. The decay temperature is therefore 



4 A flaton 4> cannot have a coupling of the form W ~ 2 i/> since it would lead to an unsuppressed quartic self- 
coupling in the potential for <f>. Such a coupling can be forbidden by an appropriate gauge symmetry. The only other 
possibility would be for cj> to be charged under some continuous gauge symmetry, which in the vacuum would be 
broken at a scale ~ vac ~ 10 10 GeV. We do not consider this possibility here. 
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T dcc ~ gjTjM^~2x 10 



/ 3/2, 



^vac 

llO 



300 MeV 



rn^ nI/ 10 10 GcV \ 
UOO/ VIOOGcV/ V ^vac / 



where 5* is the effective number of massless degrees of freedom in the universe at temperature Td cc . Note 
that parametric resonance is unlikely to be important because (p oscillates around a large vacuum expectation 
value rather than the origin. 

Renormalization, amongst other things, will split the degeneracy of Eq. ( |TT| ) amongst the various com- 
ponents of ip and ip. The renormalization will tend to make SU(3) charged components the heaviest, and 
SU(3)xSU(2)xU(l) singlet components, if they exist, the lightest Q. 

Our main unjustified assumption is now to assume that ip and ip carry opposite charge under some discrete 
(or continuous) gauge symmetry under which all other fields are neutral. This symmetry is, however, helpful 
in avoiding unwanted superpotential couplings to MSSM fields. The lightest component of ip and ip will 
then be absolutely stable and so potentially a dark matter candidate, fj For example, we could have a Zg 
discrete gauge symmetry, under which <f>, ip, and ip have charges 2, —1, and —1, respectively. This would 



guarantee Eqs. (|4|) and (10), and after cp acquires its vacuum expectation value, the Z% will be broken down 
to a Zi under which only ip and ip are charged. A simple choice of representations that satisfies the discrete 
anomaly cancellation conditions jlj] is ip = 16 + 1 and ip = 16 + 1 of SO(10). 



B. Abundance 

During thermal inflation, ip and ip will be relativistic and in thermal equilibrium. Their number density 
will therefore be 

»(T) - (14, 

where £(3) ~ 1.202, and is the number of internal degrees of freedom of ip and %p. [] If ip = 16 + 1 and 
%j) = 16 + 1 of SO(10), as in the example of the previous section, then = 136. Through the coupling of 
Eq. (|l0|), ip and ip will generate the finite temperature effective potential for <p 

V T (<P) = V ti + (|| |A^| 2 T 2 - ml) H 2 + . . . (15) 



5 Throughout this paper, we use the same notation ip for the complete representation and its lightest component 
(the dark matter). 

6 We have assumed, as is appropriate for a supersymmetric theory, that there are equal numbers of bosonic and 
fermionic degrees of freedom (7/8 = (1 + 3/4)/2). 
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Thermal inflation will therefore end at the temperature 



T C = M^L (16) 



when </> begins to roll away from <f> = 0. Shortly afterwards, the abundance of ^> and ip freezes out. Meanwhile, 
</> will continue to roll towards its vacuum expectation value \(p\ = </> V ac- Once <p acquires its vacuum 
expectation value \4>\ — (pv&c, the coupling of Eq. (^) will give ip and ip masses 

M=|A^|0 vac (17) 

The freeze-out abundance of ip and ip can be estimated as follows. First, it is important to keep in mind 
that the temperature does not drop as (p rolls from to 4> vac . This is because the time-scale for the roll-over 
is ml , which is much smaller than the Hubble time H^ 1 ~ m^ 1 <p~^ c M p i ^> . Hence T — T c throughout 
the roll-over. The freeze-out occurs as ip gains mass and begins to become non-relativistic when its thermal 
abundance is given by 

n{m^) = (j^^j ex P f-^rj (18) 

where denotes the <p dependent mass of ip, = |A^</>|, which increases as cp rolls down the potential. 
In other words, in contrast with the usual freeze-out calculation, it is that is changing with time rather 
than the temperature. The freeze-out abundance is determined by equating the annihilation rate with 
the inverse-time-scale of the problem at hand, i.e. to^. 

r^( m V>) ~ n { m i>) ( a l u l)("V) ~ m 4> ( 19 ) 

Using the fact that (<r |t>|) < 1/rrv^, and that ~ |A^| T c , it is not hard to see that the freeze-out occurs 
when m^p ~ T c , and that the freeze-out abundance is given by Eq. ( fTi] ) with T = T c , suppressed by at most 
a factor |Aw,|. So for |A^| ~ 1, and to within an order of magnitude, there is no significant net annihilation 
of ip and tp after the beginning of the roll-over, and the freeze-out occurs well before <f> reaches (p va , c . 

After the freeze-out, cp will continue to roll towards its vacuum expectation value, increasing the mass 
of ip to a large value. It can be checked that the annihilation rate by the time (p reaches the minimum is 
negligible Q 



\ |3/2 1/2* 

A | 




« Hti (20) 



7 If |A^| is small, the fifth power of \X^\ in this formula could make significant which means that there would be 
some annihilation of ip and ip. However, in this case the cosmological abundance will also be boosted up by a higher 
T c (see Eqs. @ and @). 
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Subsequently, <f> will oscillate around 4> vac . [] 

One might worry that during the oscillation, a significant amount of annihilation or production of ip and 
ip might occur if </> returns to small values. The rapid build up of gradient energy will prevent (p from 
returning to small values except in a few isolated places. In addition, <p would eventually be prevented from 
returning to small values by the Hubble expansion. Strings with walls attached are also formed, which will 
likely disappear quickly |H| . Their direct radiation into the heavy ip particles is heavily suppressed because 
^> 77i0. On the other hand, the ip particles are light in the cores of the strings, and could be created and 
trapped there. If a string loop carries a net ^-charge, it will be released in the form of (heavy) ip particles 
when the loop annihilates. If we assume each string produces of the order of one ip particle, our rough 
estimate of the ip abundance should still be valid. 

The energy density in ip and ip will then scale with that of the oscillating flaton until the flaton finally 
decays, leaving a radiation dominated universe at a temperature Td cc - The energy density of ip and ip will 
then scale with the entropy density of the universe, s. The current value of the entropy density is 

s = 2.2 x 1CT 38 GeV 3 (21) 
Finally, we wish to compare the current energy density of ip and ip with the critical density 

3# 2 = 3x 1(T 47 GeV 4 (22) 
For ip and ip to be a viable dark matter candidate we require 



|^~0.3 (23) 



Putting everything together we get 

' Pdec\ { _S0_\ f _j_ 

V ti j \s dc J \3H$ 



^ = n(T c )M lOg. * ] (24) 



where p^ cc is the energy density of the universe at the end of the flaton decay. Therefore, using p^cc = 



^dec^deo 



4 



-2 (9±\% ( m <t> \% ( 5 x 10 



W-i^M^l \ m J VlOOGeW \g^ m fM^ ) 

, / 87T 2 n(r c ) \ ( y/g^ \^\T C V / M \ ( 2m^ ac \ / m0 jVf P1 \ 
X UC(3)^T C 3J ^ 4V 2 m0 J (jA^J^ 3V ti ) [VSIX^^J [ ) 

where all factors in brackets are of order 1. We have displayed explicitly the assumed relations between the 
various quantities, e.g. M = |Ai/,|</> vac , etc. 



Note that the backreaction of the finite density of ip particles on <^'s potential is negligible. 
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III. OBSERVATIONAL CONSTRAINTS AND DETECTION 



The above simple model leaves open the question of what kind of interaction tp has with ordinary matter. 
The same is true of other production mechanisms of superheavy dark matter [Q,^). Let us go through the 
different possibilities one by one. 

Electromagnetically charged. These have been referred to as CHAMPs in the literature p"7| . At late times, 
they primarily take the form of p + ip~ (hydrogen with a heavy "electron" which has very low cross-section 
with other atoms), ip + e~ (heavy hydrogen) or ■0~He ++ e _ (ip~ bound to the helium nucleus to make another 
kind of heavy hydrogen). Various constraints exist on such particles, ranging from the absence of heavy- 
hydrogen-like atoms in water to nondetection in 7-ray and cosmic-ray detectors flTj]. By far, the strongest 
constraint appears to come from the existence of old neutron stars where only M > 10 16 GeV is allowed. 
Otherwise, a sufficient net number of r tp + particles collects in the neutron star, forms a black hole in the center 
and eats up the star on a short time scale; this is in part because the hydrogen-heavy-hydrogen scattering 
cross-section is high, given by the square of the Bohr radius in the low velocity limit a ~ 10~ 17 cm~ 2 . 
However, if the abundance of and that of tp~ bound to helium are the same in the halo, the constraint is 
weakened toM> 10 10 GeV. The conventional bound of M < 10 5 GeV would then be fatal to the existence of 
such particles. The kind of production mechanism like the one proposed here, or elsewhere, which evades the 
conventional bound, could resurrect the intriguing idea that the dark matter can be charged. But as we can 
see, significant astrophysical constraints already exist. It should be noted that the near Planck mass relic that 
will be discussed in §IV satisfies even the demanding bound of M > 10 16 GeV. The economic importance of 
such stable massive electromagnetically charged particles cannot be overestimated. For example, they could 
be used to catalyze nuclear fusion (l7) . 

Strongly charged. These have been referred to as SIMPs in the literature 2(|. Significant bounds on 
their masses, if they have significant cosmological abundances, come from nuc eosynthesis as well as the 
absence of anomalously heavy isotopes of familiar nuclei [|TJ . A systematic study of constraints from direct 
detection and the existence of old neutron stars and the Earth was made in Ref. [l9| ]. Assuming there is no 
ip-ip asymmetry, the neutron-star argument and underground plus balloon experiments provide competitive 
bounds: only M > 10 8 — 10 10 GeV is allowed for a ■i/'-proton cross-section of a ~ 1CP 30 — 10~ 25 cm -2 . As in 
the case of electromagnetically charged dark matter, the possibility of producing them without overdosing 
the universe gives such dark matter candidates a new life. 

Weakly charged. Naturally, no significant constraints exist if ip has only weak-scale interactions like the 
neutralino (i.e. a ~ 10~ 44 (to„/ GeV) 4 cm~ 2 in the large M limit, where m„ is the mass of the relevant 
nucleon). Direct detection appears rather difficult simply because the halo number density scales as M _1 , 
and the neutralino with mass ~ 100 GeV is already difficult to detect. Note that a large mass does not increase 
significantly the nuclear recoil: AE ~ M 2 m n v 2 / (M + m„) 2 where v ~ 200 km s _1 is the average halo velocity 
of these particles; in the large M limit, AE is asymptotically M- independent. Indirect detection might seem 
even more hopeless. Not only does the halo number density drop by a factor of M, the annihilation (which 
gives rise to neutrinos) rate is suppressed by M 2 according to the unitary bound. However, three opposing 
factors help us here. First, on a sufficiently long time scale, the neutrino flux from ip-ip annihilation in the core 
of the Sun or the Earth is determined not by the annihilation rate, but by the capture rate, which depends 
on the scattering cross-section with nucleons and is not heavily suppressed. Second, indirect detection works 
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by observing muons that result from the interaction of the neutrinos with the rocks of the Earth. The cross- 
section for producing muons and the range of the muons both scale up with energy, and hence the mass of 
ip. A detailed calculation taking into account these effects as well as other relevant ones will be presented in 
a forthcoming paper. 

Completely neutral. ip could be a singlet under SU(3)xSU(2)xU(l). It is of course virtually impossible to 
detect such particles, except by their gravitational effects. 

Lastly, superheavy dark matter has been postulated as the origin of ultra-high-energy cosmic-ray events at 
energies > 10 10 GeV, above the Greisen-Zatsepin-Kuzmin cut-off |22 23.1^|. The parameters of the model 



presented in the last section are sufficiently flexible to allow a mass of M ~ 10 11 GeV to explain such events. 
The ip particles cannot by themselves be the primaries because of the large mass, even if they have hadronic 
interactions |24|]2lf |. The simplest way is to have them decay into hadrons, e.g. protons, which reach the 
Earth's atmosphere. But then, one has to invoke special reasons to explain why they are not stable but 
sufficiently long-lived. 

IV. DISCUSSION 

We have shown that a simple and well-motivated model of thermal inflation naturally produces dark 
matter particles ip and ip of mass M ~ 10 10 GeV with a cosmological abundance in the correct range. As our 
mechanism is implemented by thermal inflation, which solves the moduli problem by late entropy production, 
it is robust against such dilution. 

The same general mechanism could also be applied to what one might call 'moduli thermal inflation' to 
produce near Planck mass particles at low energy scales. Moduli thermal inflation is a limit of thermal 
inflation in which the flaton is replaced by a modulus, so that roughly speaking |A^| ~ m^/Mpi in Eq. (j^). 
0vac will then be of the order of the Planck scale. In the context of string theory, it is then very reasonable 
to assume that the vacuum expectation value of the modulus, cp = (p vac , corresponds to another 'origin' in 
field space where new fields become light; for example one could have superpotential couplings of the form 
W = X x (4> — V ac) XX- Such a 'coupled' modulus would appear like an ordinary scalar field (e.g. a squark 
or slepton field) in the true vacuum. ^| The decay temperature would then no longer scale as in Eq. ( |l3] ) 
but instead could be as high as V^/ 4 . Putting these modifications into Eq. ([25| ) would also give us a value 
of f2^v, in the neighborhood of 1. However, another consequence of having </> vac ~ Mpi is that one could get 
a significant number of e-folds of non-slow- roll inflation as <f> rolls from <f> ~ to (p ~ (p vac - To avoid this 

1/2 

inflation diluting the ip particles too much, one would require > lOV^j /-Mpi and so tp vac < Mpi/10. 
This would limit the final mass of the ip particles to be M < few x 10 17 GeV. However, this is still above 
even the stringent limit on electromagnetically charged dark matter obtained in Ref. |f^] . 

Note that because moduli thermal inflation occurs at too high an energy scale to solve the moduli problem, 



It would be an excellent candidate for an Affleck-Dine field 
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this scenario would only be viable if there were no moduli problem ^| because otherwise the ifi particles would 
be diluted by another epoch of thermal inflation, or some other late entropy production, that would be needed 
to dilute the decoupled moduli produced at the end of the moduli thermal inflation. 

A related scenario could emerge from some of the more plausible models of inflation [^6| . Here ip or ip is 
the inflaton, which holds 4> at zero by the hybrid inflation mechanism [ p7| rather than thermal effects. One 
could then get dark matter in the form of charged, near Planck mass, inflatons! 
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